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MAGNETIC  -FIELD-INDUCED  ENHANCEMENT 
OF  REIATIVTSTIC  -EIECTRON-BEAM  ENERGY  DEPOSITION 


An  important  area  of  research  associated  with  the  goal  of 
relativistic-electron-beam-initiated  fusion  in  pellet  targets  has  been 
the  investigation  of  techniques  to  reduce  the  energy-deposition  length 
of  the  beam  in  dense  matter.1  Attention  has  been  devoted  to  this 
question  because  the  long  deposition  lengths  of  MeV -energy  electrons 
require  beam  powers  for  fusion2  which  axe  technologically  difficult  to 
achieve.  A reduction  in  deposition  length  in  the  target  plasma  might 
be  affected  by  the  presence  of  a strong  self-  or  external -magnetic 
field  transverse  to  the  beam-current  direction.3 

Here,  a perturbation  procedure4  is  used  to  solve  the  relativistic 
Boltzmann  equation  with  a Fokker-Planck  collision  term5  in  order  to 
characterize  the  collision&l  transport  of  relativistic  electron  beams 
in  magnetized,  high -atomic -number  plasmas.  The  analysis  is  valid  for 
arbitrary  magnetic -field  and  plasma -parameter  variations  provided  that 
macroscopic  scale  lengths  are  large  compared  to  either  a beam-electron 
scattering  length  or  gyroradlus.  The  Drdering  of  smallness  parameters 
results  in  a diffusion-like  equation  for  the  energy  distribution  func- 
tion of  beam  electrons  involving  the  electric  field  and  dynamic  fric- 
tion. Explicit  one -dimensional  solutions  are  then  used  to  determine 

Not*:  Manuscript  submitted  January  4,  1077. 


1 


the  manner  in  which  the  deposition  of  electron  energy  changes  with 
transverse  magnetic  field  strength.  Finally,  the  applicability  of  this 
mechanism  to  current  experiments  and  fusion  systems  is  discussed. 

The  equation  describing  the  momentum  distribution  function  of 
relativistic  electrons  interacting  with  a high -atomic -number  plasma 
may  be  written5 

ft + ^*^-e(s+^)‘V  = VCvs(p  )(pS?-roWp/]  +y[vE(p)p/] 

(1) 

where  Y2  = l+p2/(mc)2.  The  quantities  vg  and  vg  are  scattering  and 
energy -loss  frequencies: 

vg  = Cy/(Y2-1)3/2,  Vg.  = €Yvs  , (2) 

where  0 = 2Trn.rfc(Z2+  Z)lnA  and  € = 2/(Z  + l)  «1.  Here,  n.  is  the 
plasma  ion  density,  rQ  is  the  classical  electron  radius,  c is  the  veloc- 
ity of  light,  Z is  the  plasma  atomic  number,  and  InA  is  of  order  10. 
Since  electron  diffusive-relaxation  times  through  the  dense  plasma  are 
much  shorter  than  the  time  over  which  macroscopic  beam  parameters  vary, 
the  time-derivative  term  may  be  neglected.  In  keeping  with  the  diffu- 
sion ordering,  the  7/  and  electric -field  terms  are  considered  to  be  of 
order  €^v0.  The  magnetic -field  term  is  taken  as  comparable  to  v for 
generality. 

With  the  ordering  Just  discussed,  the  zero-order  equation  takes 
the  form 

- P x ® * Vo  = vsV  t(Pa* -PP>  * Vo3  (5) 


2 


where  / = fQ  + f + / + ...  , f ^ is  of  order  €^/Q,  fs  of  order  €/q, 
etc.  Equation  (3)  is  satisfied  only  by  a function  which  is  isotropic  in 
momentum  space,  i.e.  /o(p,x)  = fQ( P,x).  To  see  this,  p is  written  in 
terms  of  the  spherical  coordinates  (p,9,cp)  with  the  polar  axis  parallel 
to  B.  Expanding  / in  terms  of  the  surface  harmonics  Y*?(8,«p)  and 
substituting  to  Eq.  (3)  results  in  a solution  only  for  m = 1 = 0. 

Thus,  fQ  is  isotropic.  Isotropic  / is  confirmed  by  sophisticated 
particle  codes  which  have  demonstrated  that  tightly-pinched  electron 
beams  have  this  character.6 

The  first-order  equation  then  is 


=•  ( 1 eE  *'o\ 

p 7f0  Y Sf) 


= * [(p^7  -pp)  -7  M 


n 

T p 


P x b • 7p/ i 


(M 

where  b = b/b  and  = eB/m.  The  non-re lativistic  form  of  Eq.  (4)  has 
been  solved.4  The  solution  follows  from  inspection  of  reference  4. 


where 


(5) 


(6) 


and  a = no/2Vvg.  The  equation  for  /o(p,x)  is  determined  by  integrating 
Eq.  (l)  over  the  momentum-space  solid  angle  and  keeping  terms  to  lowest 
significant  order.4  This  procedure  results  in 
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(7) 


In  the  limit  of  Qq  = 0,  Eqs.  (5)  and  (7)  reduce  to  published  results.7 
The  interesting  macroscopic  quantities  derivable  from  solution  of 

Eqs.  (5)  and  (7)  are  the  particle  flux  ? and  heat  flux  q.  These  are 


defined  by 


* 

( 

l 

q 

J 1 

mc2(Y-l) 

(8) 


Taking  the  divergence  of  Eq.  (8),  substituting  from  Eq.  (7)>  and 


integrating  by  parts  over  p results  in  the  relations 


7 . ? = - 4TT€ns(mc)3/o(o,x)  (9) 

Q = -7»q  = eE*T  + 4tt  € ng(mc)3*  me2  jY2/QdY  . (10) 

The  particle  "sink"  in  Eq.  (9)  represents  the  merging  of  beam  elec- 
trons with  the  plasma  background  once  they  have  been  slowed  by  dynamic 
friction  to  low  energies.  The  quantity  Q is  the  volumetric  heating 
rate  of  plasma  due  to  both  ohmic  and  collisional  beam-energy  losses. 

One -dimensional  solutions  in  a slab  geometry  from  which  ? and 
Q can  be  determined  explicitly  are  now  considered.  A monoenerget ic 
(Y  = Y ) electron  beam  of  particle  flux  iQ  is  assumed  incident  on  a 
uniform  plasma  occupying  the  half- space  x * 0.  A uniform  magnetic 
field  of  arbitrary  strength  aligned  parallel  to  the  plane  x = 0 is 
imbedded  in  the  plasma. 


k 


Equation  (9)  in  one  dimension  states  that  the  beam  flux  decays 

with  x so  that  in  order  to  maintain  charge  conservation,  a plasma- 

electron  current -density  j(x)  = e($-$Q)  must  flow.  An  electric  field 

E(x)  = T| J , T|  being  the  plasma  resistivity,  must  then  exist.  The 

importance  of  E to  evaluation  of  fQ  is  determined  by  comparing  the 

magnitude  of  the  7/q  and  electric -field  terms  in  Eq.  (4).  it  is  now 

shown  that  E is  negligible  for  problems  of  interest,  i.e.  that 

eE  « in^c^y^x  1>  where  x is  the  scale  over  which  / varies.  Since 

E eTl$0>  and  X ^ c/ for  electron  energies  of  interest®,  the  above 

condition  is  satisfied  if  eT)J  « €^n  c2  Y n_/e.  Substituting  numerical 

o o o S 

values  and  assuming  classical  resistivity  it  is  found  that  the  electric- 
field  is  negligible  when 


e$  « 3 x 10: 
o 


Amps/ cm* 


(H) 


In  this  expression,  n^/n,,  is  the  ratio  of  plasma-ion  to  solid  density, 

9 is  the  plasma  temperature  in  eV  and  is  it’s  ionization  level.  The 
numerical  coefficient  is  appropriate  for  a gold  plasma.  This  condition 
is  satisfied  for  electron  beams  of  interest  so  that  the  electric  field 
can  be  ignored  for  determination  of  fQ. 

With  the  simplifications  discussed,  Eq.  (7)  reduces  to 


*'  + [<^f + 5?]  &**'.>  ■ 0 M 

where  a?  = c2/6€fl2.  This  relation  takes  the  form  of  a simple  diffusion 
equation  with  the  change  of  variable 
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t(Y)  = 


dY 


(13) 


for  which  the  second  term  in  Eq.  (12)  reduces  to  -9(y2/0)/9t. 
The  solution  corresponding  to  monoenergetic  /q  at  x = 0 is8 


x 


with  determined  from  the  solution  of  Eq.  (9)  and  the  condition  that 

*(o)  = *Q. 

In  the  limit  of  large  t can  be  approximated  by 

t “•  4r|(Yo.v)/^  for  which  Eq.  (10)  can  be  integrated  to  give  the 
deposition  profile  in  a strong  magnetic  field 


1 = TT^k*erfc(kx)  ; 

J Qdx 


0 

o 

1*ns(vo'1)^ 


(14) 


o 

3/ 

This  form  is  accurate  for  fio/fig  > 8y^/(y^-1) '2.  For  lower  magnetic 
fields, Q(x)  is  determined  by  numerical  integrations  of  Eqs.  (10)  and 
(13).  Results  are  displayed  in  Fig.  1 for  YQ  = 3.  The  case  = 3 

results  in  a profile  indistinguishable  from  Eq.  (l4) . 

The  magnetic -field-induced  enhancement  in  peak  heating  rate  (or 
reduction  in  deposition  length)  above  that  due  to  dynamic  friction 
alone  can  be  estimated  from  Eq.  (12)  by  identifying  the  character- 
istic deposition  length  A 


A(no)  = a 


Yo 

(Y§-D3 


(15) 
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Then 


tf^y2-!)3 

ov  o ' 

4nfv4 

S o 


(16) 


The  numerical  results  show  that  this  relation  is  valid  to  within  about 

10$  for  all  values  of  Q for  1 MeV  incident  electrons. 

o 

It  is  of  interest  to  determine  whether  magnetic  enhancement  of 
energy  deposition  is  operational  in  current  experiments  and  scalable  to 

fusion  conditions.  This  is  best  accomplished  by  numerically  solving  Eq. 

, . =t 

(7)  with  M determined  from  the  azimuthal  self- field  of  the  pinched  beam 
or  by  employing  Monte -Carlo  techniques.3  However,  when  the  beam  radius 
R exceeds  A,  Eq.  (l6)  with  determined  by  B„  can  provide  an  estimate. 
Magnetic  enhancement  is  important  when  A > 1.  This  condition  is  satis- 
fied in  the  target -plasma  blow-off  when  Bg  > 2 X 10*^/ ng  gauss  for  a 
1 MeV  beam  incident  on  gold.  In  this  case,  R exceeds  A when  the  beam 
current  exceeds  50  kA.  The  condition  A > 1 is  then  satisfied  in  current 
experiments9 >10  provided  that  the  plasma  is  sufficiently  resistive  to 
allow  B.  to  diffuse  a distance  A into  it.  Thus,  a necessary  ccndition 
for  deposition  enhancement  due  to  self -fields  is  t > td  = 4ttA2/T]c2, 
where  t is  the  beam  duration.  For  A > 1 in  gold,  this  condition  may  be 

3/ 

written  t > 1280/2/z^B2.  Evaluating  rD  for  9 = 10  eV,  Z^  = 5>  and 
Bg  = 2 x 10s  g results  in  a value  of  20  ns.  Note  that  although  the 
diffusion -time  condition  is  satisfied  in  current  experiments,  it  is  not 
sufficient  to  insure  B-field  penetration  since  plasma  expansion  can 
convect  field  lines  away  from  the  dense  region  where  enhanced  deposi- 
tion is  desired.  Simultaneous  solution  of  the  plasma  MHD,  deposition, 
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and  magnetic -diffusion  equations  in  a realistic  geometry  is  required  for 
a definitive  answer  to  the  role  of  magnetic  enhancement  in  current 
experiments.  However,  at  the  elevated  temperatures  required  for  fusion- 
pellet  targets,  these  considerations  indicate  that  it  is  unlikely  that 
sufficient  self-fields  would  penetrate  the  plasma  unless  anomolous 
resistivity  were  present  or  unless  a significant  fraction  of  the  field 
could  be  made  to  penetrate  before  plasma  heating  occurred.  This  last 
situation  might  occur  if  the  beam  were  focussed  onto  the  pellet  only 
after  its  current  has  risen  close  to  the  maximum  value,  as  is  done  in 
certain  experimental  configurations.®  In  that  case,  the  unpenetrated 
portion  of  the  field  must  confine  the  heated  plasma  in  a way  which  does 
not  disrupt  synmetric  implosion  of  the  pellet. 
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